arXiv:1509.03694vl [astro-ph.GA] 12 Sep 2015 


Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 15 September 2015 


(MN Wr^ style file v2.2) 


Variability of broad absorption lines in a QSO SDSS 
J022844.09+000217.0 on multi-year timescales 

Zhi-Cheng He, Wei-Hao Bian*, Xiao-Lei Jiang & Yue-Feng Wang 

Department of Physics and Institute of Theoretical Physics, Nanjing Normal University, Nanjing 210023, China 


15 September 2015 


ABSTRACT 

The variability of the broad absorption lines is investigated for a broad absorption line 
(BAL) QSO, SDSS J022844.09+000217.0 (z = 2.719), with 18 SDSS/BOSS spectra 
covering 4128 days in the observed frame. With the ratio of the rms spectrum to 
the mean spectrum, the relative flux change of the BAL-trough is larger than that 
of the emission lines and the continuum. Fitting the power-law continuum and the 
emission line profiles of C iv A1549 and Si IV A1399, we calculate the equivalent width 
(EW) for different epochs, as well as the continuum luminosity and the spectral index. 
It is found that there is a strong correlation between the BAL-trough EW and the 
spectral index, and a weak negative correlation between the BAL-trough EW and the 
continuum luminosity. The strong correlation between the BAL-trough EW and the 
spectral index for this one QSO suggests that dust is intrinsic to outflows. The weak 
correlation between the BAL variability and the continuum luminosity for this one 
QSO implies that the BAL-trough variation is not dominated by photoionization. 

Key -words: galaxies:active—galaxies:nuclei—quasars:absorption lines 


1 INTRODUCTION 


Broad absorption line quasars (BAL QSOs) exhibit broad, 
ultraviolet (UV) line absorption troughs. The traditional 
BAL QSOs were quantified by balnicity index (BI) with 
velocity width larger than 2 000 km s~^ and outflow veloc¬ 
ity of 3000-25000 km s”’^ dWevmanrJ Il99ll l. There are a 
large number of broad absorption features within 3000 km 
s“^ and beyond 25000 km s“^ , as well as features with 
widths less than 2000 km s~^ , and newer ways of quan¬ 
tifying the in dex (e.g. iHall et ^ I 2 OO 2 I : I Trump et al.l I 2 OO 6 I : 
iGibson et al.l l2009ll. BAL troughs are present in about 10- 
40% of QSOs (e.g. Hewett &c Foltz gOOsI : FlVump et aDl2006l : 
iGangulv et al.ll2()07 : Allen et al.ll2011 1. produced by absorp¬ 
tion from high-ionization lines such as Si iv 1399, G iv 
1549, C III] 1909 (known as HiBAL), and low-ionization lines 
such as Mg II2799 (known as LoBAL). The fraction of BAL 
QSOs strongly depends on the dehnition of BAL QSOs (e.g. 
iHewett fc Foit j[2003l : iTrump et al.ll2005l , so drawing physi¬ 
cal conclusions from population statistics must be done with 


BAL troughs are tho ught to be th e strongest observed 
signatures of QSO winds (lFabianll2012h . The observation of 
BALs is commonly thought to be the result of passing a disk 
wind along the line of sight (e.g. [Murray et al.l[l995l : lElvisI 
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I 2 OOOII . If the wind is on the line of sight, a QSO appears as a 
BAL QSO. That detection of a wind should be orientation- 
dependent is very similar to the case of other structures 
in QSOs, such as the broad line region (BLR) or jet (e.g. 
lUrrv fc Padovanlll995f) . 


However, it was found that there is no signihcant 
difference between radio spectral index distributions for 
radio-loud BALs and radio-loud non-BAL s, indicating 
they h ave similar ranges of viewing a n gles (iBecker et al.l 
2 OOOI: Montenegro-Montes et al.] l2008l : iFine et al.l I 201 H : 


Bruni et al.l 2012ll . No correlations exist between outflow 


properties and the orientation, suggesting that BAL winds 
along any line of sight are driven by the same mechanisms 
(DiPompeo et al. 2012). There are also B AL QSOs that are 
seen very nearly along the radio jet axis llZhou et al.ll20oJ : 
iGhosh fc Punsl'Jl2007l) . Another explanation f or BALs is an 
orien t ation-independent evolu tion effect (e.g. I Gibson et al.l 
l2008l : IZubovas fc Kin j l2013ll . BALs outflow is possibly 
caused by the expulsion of gas and dust by galaxy colli¬ 
sion as an evolutionary stage of AGN (e.g. IVoit et m]I 1993I : 
iGregg et al.ll2006l : iGibson et al.ll2008h . 


The disk wind is believ ed to come from t he central 
accretion disk in QSOs (e.g. [Murray et al.l[l9^ . The de¬ 
pendence of wind properties (such as the trough EW, the 
maximum velocity, the fraction of BAL QSOs) on the phys¬ 
ical properties of QSOs (such as the luminosity, the He ll 
A1640 A EW, and the UV continuum slope) has been inves- 
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tigated by man y authors, and they found some correlations 
between them dLaor fc Brandd I 2 OO 2 I : iGangulv et alJ l2007l : 
I Baskin et alJ 120131 '). I~ was found that the maximum out¬ 
flow velocities increase with both the bolometric luminosity 
and the blueness of the spectral slop e, suggesting the idea 
of radiation pressure -driven outflows (iLaor 8z Brandtl|2002l : 
IGangulv et al.l[2007^ . 

BAL troughs often vary in equivalent width (EW) 
and/ or shape over rest-f r ame timescales of months to years 
(e.g. I ^rlow et H] Il992l : iGibson et al.lI 2 OO 8 I : iFiliz Ak et al.l 
I 2 OII I 2 OI 3 I) . It was found that the fractional EW change 
increas es with rest-frame timescale over the ra nge 0.05 - 5 
yr fe.g. iGibson et al]l2008l : IPiliz Ak et al.ll201^') . The BAL- 
trough variation is believed to be driven by changes in 
velocity structure, c overing factor, or ionization level (e.g. 
iFiliz Ak et al.ll2012ll . It was suggested that BAL variation 
arises from chang es in the amount of ’’shielding gas” along 
the line-of-sight llFiliz Ak et al.l l2013l L. Such changes are 
likely according to disk-w ind simulations fe.g. IProga et al] 
l200d : I Murray et al.l 1 19951 1 . These changes coul d ultimately 
be re lated to the central accretion disk (e.g. IProga et al.l 
l2000h . Changes in the column density of shielding gas (BAL 
trough EW variability) can relate to the level of ionizing 
luminosity reaching the BAL wind. 

Two kinds of investigations of BAL-trough variabil¬ 
ity have typically been undertaken. One is to use two- 
epoch spectra in a BAL sample; the ot her is to use multi¬ 
epoch spectra for one single object (e.g. Gibson et_ajJ^l2008l : 
IFiliz Ak et al.l[2013l : ICaDelluDO et al.ll201ll . 2012I . 2013l l. Us¬ 
ing a two-epoch sample of 13 BAL QSOs overlapping be¬ 
tween the Very Large Array (VLA) Faint Images of the 
Radio Sky at Twenty-Centimete r (FIRST) Su r vey a nd the 
Sloan Digital Sky Survey fSDSSl. [Gibson et al.l (I 2 OO 8 I I found 
that BALs tend to vary on multi-year timescales in veloc¬ 
ity regions that are a few thousand kilometers per second 
wide. They also found that BAL-trough variations do not 
appear to be correlated with variations in the observable 
continuum, implying that the BAL-trough va riation is not 
dominated by photoionizat io n changes (e.g. Barlow et al.l 


198E: l^pellupo et al.l[201^ . ICapellupo et"^! ( 201ll l2012l . 
201311 did extensive variability studies on a small BAL QSO 


sample with longer timescales and multiple epochs. They 
found that the variability occurs typically in only portions 
of the BAL troughs; the components at higher outflow ve¬ 
locities are more likely to vary than those at lower veloci¬ 
ties, and weaker BALs are more likely to vary than stronger 
BALs. The shortest timescales constrain the location of the 
outflo wing gas. With a larger SDSS sample, IFiliz Ak et al.l 
1 I 2 OI 3 II found that there is a relation between the luminosity 
and EW variability on moderate timescales, although they 
did not regard this case as strong evidence for luminosity 
dependence. 

The BAL-trough variability for one single object is very 
important in understanding of the origin of the BAL QSOs. 


various data releases (e.e. Trumo et al. 

I 2 OO 6 I: IGibson et al.l 

l2009l:IShen et al.ll201ll:IParis et al.ll2013l 

). We selected a sin- 


gle QSO, SDSS J022844.09+000217.0, which has the high¬ 
est number of observations in the time domain, to investi¬ 
gate the relation between the wind and QSO properties. §2 
presents the data. §3 gives the data analysis. §4 contains 
our results and discussion. The Summery is given in the 
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Figure 1. A typical single spectrum (top panel), the average 
spectrum (second panel), the root-mean-square (rms) spectrum 
(third panel), the ratio of the rms spectrum to the average spec¬ 
trum (i.e., the relative flux change, bottom panel) from total 18 
spectra for SDSS J022844.09+000217.0. The dotted lines show 
the wavelengthes of Si iv A 1400 A and C IV A 1549 A. 


last section. Throughout this work we use a cosmology with 
Ho = 70km s“^Mpc“^, D,m = 0.3, and Da = 0.7. 


2 SDSS/BOSS SPECTRAL DATA 

The Baryon Oscillation Spectroscopic Survey (BOSS), part 
of SDSS-III, uses the dedicated 2.5-m wide-field telescope at 
Apache Point Observatory near Sacramento Peak in South¬ 
ern New Mexico to conduct an imaging and spectroscopic 
survey for about 1.5 million lumin ous galaxies as well a s 
about 160,000 quasars at 2 > 2.2 llEisenstein et al.]l2fllll ). 
With respect to the original SDSS, the wavelength coverage 
of BOSS spectra changes from 3800 A -9200 A to 3600 A 
-10400 A, as well as the fiber numbers per plate from 640 to 
1000. And the BOSS uses new fibers (1000 rather than 640 
per plate), with small er holes (2” rathe r than 3”). 

With SDSS DR7, IShen et ^ (I 2 OIII ) gave a compilation 
of properties of the 105783 QSOs. There are 6214 BA L QSO s 
in their SDSS DR7 sample (see also iGibson et akl l2009l ). 
There are 456 G iv BAL QSOs with at least two-epoch spec¬ 
troscopic observations. We searched for observations of these 
456 C IV BAL QSOs in SDSS DRIO and found more epochs 
of spectroscopy for each object. The one BAL QSO with the 
most observations, SDSS J022844.09+000217.0 (z = 2.719), 
has 18 high signal-to-noise ratio (S/N) SDSS/BOSS spectra, 
covering 4128 days (about 11 years) in the observed frame 
(Table 1), i.e., about 3 years in the rest frame. For SDSS 
spectra, the S/N in the r-band is about 17-20. The BOSS 
spectra have longer exposure times than the SDSS sample, 
and have a better S/N about 22-36. In Table 1, the informa¬ 
tion for these 18 spectra is listed, such as the survey name, 
MJD, Plate, Fiber ID, S/N. 
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Table 1. The properties of 18-epoch spectra for SDSS J022844.09+000217.0. Col(6-9): in units of A. Col(lO): measured from ISOOA to 
2400A,in units of lO'^^erg s“^A“^. Col(12): in units of km s“^ . EWi is measured based on a power-law continuum. EW2 is measured 
based on a pseudo-continuum of power-law continuum plus emission lines. 


Survey 

MJD 

Plate 

Fiber 

SN(r) 

£;W'i(C IV ) 

EW2{C IV ) 

EWi{Si IV ) 

EW2{Si IV ) 

Leant 

a 

^max 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

(10) 

(11) 

(12) 

SDSS 

51817 

406 

35 

17.1 

12.86 ±0.79 

18.14 ± 1.23 

8.59 ±0.55 

9.48 ± 0.62 

3.69 ±0.16 

-0.76 ±0.02 

12738 

SDSS 

51869 

406 

37 

19.6 

15.33 ±0.91 

19.32 ± 1.12 

9.77 ±0.49 

10.66 ±0.54 

3.20 ±0.13 

-1.16 ±0.03 

13124 

SDSS 

51876 

406 

37 

17.0 

14.35 ± 1.01 

18.13 ± 1.24 

8.90 ±0.58 

10.26 ± 0.67 

3.50 ±0.15 

-0.79 ±0.03 

12738 

SDSS 

51900 

406 

36 

17.0 

14.96 ±0.80 

19.58 ± 1.32 

9.29 ±0.57 

10.08 ± 0.63 

3.72 ±0.17 

-0.74 ±0.03 

12738 

SDSS 

52200 

705 

431 

15.6 

15.62 ± 1.21 

20.07 ± 1.52 

9.23 ±0.63 

10.29 ±0.70 

2.99 ±0.16 

-0.75 ±0.03 

13317 

SDSS 

52205 

704 

640 

16.6 

16.00 ± 1.03 

19.57 ± 1.35 

8.93 ±0.52 

9.64 ±0.57 

2.94 ±0.16 

-0.90 ±0.03 

13124 

SDSS 

52238 

406 

39 

19.3 

16.70 ± 1.11 

20.55 ± 1.35 

9.20 ±0.62 

9.82 ± 0.67 

3.59 ±0.13 

-0.77 ±0.03 

13703 

BOSS 

55179 

3615 

780 

36.6 

17.02 ±0.54 

22.67 ±0.71 

10.97 ±0.31 

12.32 ± 0.35 

3.03 ±0.07 

-0.32 ±0.02 

13703 

BOSS 

55181 

3647 

778 

28.6 

16.10 ±0.67 

22.33 ±0.91 

10.16 ±0.38 

11.33 ±0.43 

3.16 ±0.09 

-0.32 ±0.02 

13703 

BOSS 

55208 

3615 

780 

29.0 

16.24 ±0.59 

21.60 ±0.77 

10.18 ±0.34 

11.16 ±0.38 

2.67 ±0.07 

-0.67 ±0.02 

13703 

BOSS 

55209 

3744 

575 

22.3 

15.15 ±0.83 

20.35 ± 1.10 

10.08 ±0.51 

11.03 ±0.49 

2.90 ±0.09 

-0.62 ±0.03 

13896 

BOSS 

55241 

3647 

778 

33.8 

16.46 ±0.57 

22.56 ±0.77 

10.33 ±0.33 

11.47 ±0.35 

3.19 ±0.07 

-0.46 ±0.02 

14089 

BOSS 

55445 

3615 

770 

27.7 

15.94 ±0.73 

21.61 ±0.97 

9.97 ±0.41 

11.17 ±0.36 

2.87 ±0.09 

-0.41 ±0.02 

14089 

BOSS 

55455 

4238 

736 

30.8 

16.40 ±0.63 

22.09 ±0.83 

9.56 ±0.32 

10.50 ±0.42 

3.13 ±0.09 

-0.73 ±0.02 

14089 

BOSS 

55476 

3647 

786 

34.0 

16.58 ±0.59 

22.54 ±0.78 

10.34 ± 0.32 

11.49 ±0.40 

3.43 ±0.09 

-0.54 ±0.02 

13896 

BOSS 

55827 

3647 

780 

37.0 

17.69 ±0.67 

22.80 ±0.86 

10.86 ±0.38 

12.07 ±0.43 

3.22 ±0.10 

-0.62 ±0.02 

14089 

BOSS 

55856 

3615 

776 

32.5 

17.83 ±0.80 

22.94 ± 1.04 

11.27 ±0.44 

12.74 ±0.50 

3.16 ±0.12 

-0.40 ±0.02 

14089 

BOSS 

55945 

3647 

738 

35.4 

16.11 ±0.65 

22.14 ±0.90 

10.67 ±0.40 

11.95 ±0.45 

3.53 ±0.12 

-0.44 ±0.02 

13896 


Table 2. Summary of the Spearman correlation coefficients: EWi is measured based on a power-law continuum. EW2 is measured based 
on a pseudo-continuum of a power-law continuum plus emission lines. EW^ is measured within the high-velocity part. EWi is measured 
within the low-velocity part. The first two lines are for the C iv BAL-trough, and the last two lines are for the Si iv BAL-trough. The 
value in brackets is the probability of the null hypothesis. Leant is measured from 1300A to 2400A. 



EWi 

EW 2 

EWh 

EWi 

C IV Leant 

a 

-0.19 (0.443) 
0.45 (0.06) 

-0.21 (0.399) 
0.77 (0.00018) 

-0.22 (0.39) 
0.64 (0.0042) 

-0.28 (0.26) 

0.87 (0.000003) 

Si IV Leant 

cx. 

-0.24 (0.33) 
0.78 (0.00013) 

-0.22 (0.37) 
0.81 (0.00004) 

-0.35 (0.15) 
0.75 (0.0003) 

-0.22 (0.39) 

0.82 (0.00003) 


3 ANALYSIS 

3.1 Mean and RMS spectra 

The mean spectrum from these 18 spectra is calculated as 
follows: 

f fAl - V 

i=l 

where the sum is taken over the N = 18 spectra. From the 
mean spectrum with high S/N, the absorption region and 
emission region are more obvious, which can help us choose 
suitable continuum windows. The root-mean-square (rms) 
spectrum from the 18 spectra is determined by: 

I" iV 1/2 

/rms (A) = -j ^ ', [fiW ~ /mean (A)] I 

where the sum is taken over the N = 18 spectra (e.g., Kaspi 
et al. 1999). For the spectrum at MJD = 52205, the data 
from 1660 A - 1798 A (rest wavelength) are not used. The 
mean spectrum and the rms spectrum are shown in Fig. 1., 
as well as a single spectrum (top panel). 

From the mean spectrum (second panel in Fig. 1), the 


broad absorption troughs in Si iv and C iv are obvious, as 
well as some narrow absorption lines. From the rms spec¬ 
trum, the emission region varied significantly larger than 
the absorption region (third panel in Fig. 1). The rms-to- 
mean ratio spectrum, i.e., the relative flux change, is shown 
in the bottom panel of Fig. 1. The variance fraction from 
these 18 spectra increases as the frequency increases, from 
about 10% on the red side to about 15% on the blue side. 
The relative flux change is larger for the absorption trough 
than for the continuum and the emission lines. 

In Fig. 2, we also show all 18 spectra normalized at 2000 
A. There is clear variability of the spectral slope during this 
long timescale of about 11 years. 


3.2 Fitting the spectrum 

In order to investigate the relation between the BAL-trough 
and QSO properties, we fit the continuum and the emission 
line profile s as follows 417 et al +00^ : IZhang et al.l l2010l : 
iBian et aDl2012l : ICapellupo et al.ll2012l '): 

(1) Fitting the continuum. All the observed spectra 
are corrected for Galactic extinction using an Av value 
of 0.1575, assuming the extinction curve of Cardelli et al. 
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Rest Wovelength (A) Rest Wavelength (A) 

Figure 2. All 7 SDSS (left) and 11 BOSS (right) spectra normalized at 2000 A with a smoothing box of 20 A. The power-law continuum 
is plotted in this figure. The black horizontal lines are the continuum windows for the power-law fit. 

spec-3647-55241-0778,fits z=2.717 spec-3647-55241-0778_cln,fits 




Figure 3. An example of the power-law continuum fit (left) and the emission line fit (right). Left top: The extinction-corrected rest- 
frame spectrum is shown by the black lines. The green dots are the initial continuum windows. The blue line is the power-law continuum. 
Left bottom: the residual spectrum. Right top: The multiple Gaussian components are in blue and the sum of these is in red. The fitting 
window is in blue. The horizontal pink dots are the final data used in emission line fit. Right bottom: the residual spectrum. 


(1989; IR band; UV band) and O’Donnell (1994; optical 
band) with Rv ~ 3.1. And they are converted to the 
rest frame by redshift of z = 2.717. It is common to use 
the power-law formula, f\ oc A“ {f„ oc to ap¬ 

proximately fit the QSO continuum spectrum. Other con¬ 
tinuum models, such as a polynomial function or a red- 
den ed power-law, have a l so been used in pr evious stud¬ 
ies (iLundgren et al.l l2007l : I Gibson et al.l [200^ . We fit the 
power-law continuum iteratively in the "continuum win¬ 
dows”, which are known to be relatively free from strong 
emission lines. The fit is performed by minimizing x^- 
From Fig. 1, we select the following continuum windows: 
1290-1300, 1685-1715, 1965-2000 A in the rest frame (e.g 
ForsteretalJ 201ll : Vanden Berk et ^ 200 ll : iGibson et al.l 


20081 : Bian et ahl 20121 : Baskin et al.l 20131 ). At each itera¬ 


tion, we ignore any spectral bins that deviate by more than 
3a from the previous iteration’s continuum fit. That itera¬ 
tive approach would exclude additional spectral regions that 
contain broad emission or absorption features, especially for 
BAL QSOs. For this object, we don’t fit the UV Fe ll emis¬ 
sion. The power-law continuum is also shown in Fig. 3. An 


example of the continuum fit and the residual is shown in 
the left panel of Fig. 3. 


(2) Fitting the broad emission lines. Gonsidering the 
broad wing in the broad emission line profiles, we use three 
sets of two Gaussians to simultaneously model emission line 
profiles for Si iv A 1399, C iv A 1549, C iii] A 1909 in the 
continuum-subtracted spectra. Each emission line is fit iter¬ 
atively. At each iteration, the ’’absorbed” bins are ignored 
when they are more than 2.5cr below the previous model fit. 
In order to eliminate absorption trough impact on the fitting 
of the emission lines of the C iv , Si iv , C ill] , especially 
in the blue wings of these lines, the fitted weights of the 
blue wings for these emission lines are set to half of normal. 
And the weights of the red wings are set to twice normal. 
An example of the line fit and the residual is shown in the 
right panel of Fig. 3. The horizontal pink dots are for the 
line fitting. The red lines show the sum of the blue Gaussian 
profiles for the different lines. 
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3.3 The BAL EW and the maximum velocity of 
outflow 

Two methods are used to calculate the BAL-trough EW. 
One is based on the power-law continuum, the other is based 
on the pseudo-continuum of the power-law plus the emission 
lines. The latter includes the BAL correction from the emis¬ 
sion lines. The equivalent widths of the BAL-troughs for Si 
IV and C iv are calculated as follows: 

EW= [{1- 

J JconyA) 

The integration is done for fobviX) < fcon{X). From Fig. 5, 
for the Si iv A 1399 BAL, its EW is integrated from 1340 A 
to 1400 A, and for the C IV A 1549 BAL, its EW is integrated 
from 1450 A to 1550 A. The EW results are shown in Cols. 
6-9 in Table 1, where the subscript of 1 is for the case based 
on the power-law continuum [EWi), and the subscript of 2 
is for the case based on the pseudo-continuum of the power- 
law plus the emission lines [EW 2 )- The latter EW is larger 
than the former one. The EW 2 are about 10% larger than 
the EWi for Si iv , while the difference is about 25% for C 

IV . 

We also measure the maximum velocity of outflow. We 
searched each spectrum for absorption troughs at a level of 
10% below the continuum as given in the traditional defini¬ 
tion of a BAL. The maximum velocity corresponds to the 
shortest wavelength for a given trough with at least 10 A 
deeper than 10% below the continuum. The result is listed 
in Col. (12) in Table 1. 

For the power-law continuum, fx = / 2 ooo( A/2000A) “, 
the error for the power-law continuum at the wavelength of 
A by error propagation is: 

<5(/=on) = )" + (In A - In 2000)^^^ 

V /2000 

where the errors of S{f 2000 ) and 5a are given in the power- 
law fitting. The error for the BAL EW is measured as fol¬ 
lowed: 

siEW) = 

W ^ Jeon Jobs Jeon 

where 5{fobv) at A is the error in the SDSS/BOSS spectrum. 



Rest Time(years) 


Figure 4. The light curves for C iv BAL EW, the continuum 
(measured from 1300 A to 2400 A) and ol. The continuum is in 
units of 10^^erg 



A T (years) 


Figure 5. C IV BAL EW variation versus timescale for two- 
epoch spectra in the rest frame. The cyan line is the standard 
deviation derived from the data, calculated using a sliding window 
containing 20 time-ordered data points. The standard deviation 
of EW variation increases with increasing rest-frame timescale. 


4 RESULTS AND DISCUSSION 

4.1 Relation between the C iv emission line 
luminosity and the continuum 

There is a correlation between the C iv luminosity and the 
continuum, where the continuum is the mean power-law lu¬ 
minosity between 1300 A and 2200 A. The Spearman coeffi¬ 
cient R is 0.67 with the the probability of the null hypothesis 
of PnuU = 0.002. It is consistent with the photoionization 
scenario and classic Baldwin effect for the emission lines in 
QSOs. 


4.2 The light curves for BAL-trough C iv EW 

Fig. 4 gives the light curves of C iv EW, as well as the contin¬ 
uum, and the spectral index a (from top to bottom panel). 


The shortest time interval is one day, and the longest time 
interval is about 11 years in the observed frame. With time, 
there is an increasing trend of C iv EW and the spectral 
index, and a decreasing trend of the continuum. The largest 
EW changes are about 25% (from 18.1 A to 22.7 A) for C 
IV and 28% (from 8.6 A to 11.0 A) for Si iv . 

It was found that C iv BA L-trough variability 
is larger for longer timescales (e.g . Gibson et al.l 120081 : 
ICapellupo et al.ll2012l : lFiliz Ak et al.ll201^ ~). We also give the 
EW variation versus the time interval for two-epoch spec¬ 
tra (Figure 5). For longer time intervals, the EW variation is 
larger for th is object. That result is consistent with previous 
studies (e.g. iFiliz Ak etTll2013l l. 

We also find that there is a strong correlation between 
the C IV BAL-trough EW and the Si iv BAL-trough EW. 
The Spearman correlation coefficient is 0.7 (Pnuii = 0.003) 
for EWi and 0.8 (Pw; = 7 x 10"®) for EW2. 
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a a 


Figure 6. The continuum versus the spectral index a. The con¬ 
tinuum is measured at 1300 A (left top), at 1700 A (right top), at 
2000 A (left bottom),and from 1300 A to 2400 A (right bottom). 
The correlation becomes stronger when the continuum luminosity 
is measured for the bluer bands. 

4.3 The continuum versus the spectral index a 

For the slope variation in QSOs, some investigations found 
that spectra of QSOs with low redshift are bluer during thei r 
brighter phases (e.g. IVanden Berk et al.l 1200 ll : |Pul 20061) . 
However, with two-epoch variation, it was recently found 
that the spectra of half of the QSOs appear re dder during 
their brighter phases, especially for high-z QSOs jBian et al. 
I2OI2I) . From Fig. 1, we find that the variability is larger for 
the blue side in the rms or rms/mean spectra. The slope a 
changes about 260% (from -0.32 to -1.16, Table 1). In Fig. 6, 
we show the relation between the continuum and a. In order 
to avoid accidental conclusions, we measured four different 
continuum luminosities from the mean power-law luminos¬ 
ity at different monochromatic wavelengths or wavelength 
interval: (1) 1300 A to 2400 A; (2) 1300 A; (3) 1700 A; (4) 
2000 A. The continuum luminosity (from 1300 A to 2400 
A) changed by about 38%. The Spearman correlation coef¬ 
ficients between the continuum luminosity and a are -0.69 
(0.001), -0.36 (0.14), -0.001 (0.997), -0.30 (0.23) (from top 
to bottom and from left to right in Fig. 6), where the values 
in brackets are the probabilities of the null hypothesis. The 
correlation becomes stronger when the continuum luminos¬ 
ity is measured for the bluer bands. This correlation shows 
that this QSO has a steeper UV slope when its luminosity 
increases. This kind of negative correlatio n is consis tent with 
some previous studies for PG QSOs (e.g. |Pull2006l) . 

4.4 EW versus the spectral index a 

Fig. 7 displays the relation between the BAL-trough EW and 
the spectral index a for C iv (top) and Si iv (bottom). The 
blue squares denote the EWi data, i.e., based on the power- 
law continuum. Red circles denote the EW2 data, i.e., based 
on the pseudo-continuum of the power-law plus the emission 
lines. Due to the correction of the emission lines, the value of 
EW2 is larger than EWi. The correlation between the EW 
and a is strong. The correlation coefficient is 0.45 (0.06) 
for C IV EWi, 0.77 (0.00018) for C iv EW2, 0.78 (0.00013) 
for Si IV EWi, and 0.81 (0.00008) for Si iv EW2, where the 
values in brackets are the probabiliies of the null hypothesis. 



a 



a 


Figure 7. The BAL-trough EW versus spectral index a for the C 
iv BAL (top), and for the Si IV BAL (bottom). The blue square is 
for the case based on the power-law continuum fit. The red circle 
is for the case based on the pseudo-continuum fit of the power-law 
plus the emission lines. 


With the correction from the emission lines creating deeper 
troughs, stronger correlations are found for the C iv and 
Si IV BAL-troughs. We list these correlation coefficients in 
Table 2. 

This result shows that the BAL-trough EW becomes 
larger when the UV slope becomes flatter, implying that 
the outflow has the effect of r edden ing the observed spec¬ 
trum. Recently, I Baskin et al.l ll2013l ) used the mean SDSS 
spectra constructed by selecting on various QSO properties 
to investigate the outflow dependence on those properties. 
They found that the BAL-trough becomes deeper for QSOs 
with low He ll A1640 A EW or flatter UV slope. Our result is 
consistent with theirs. The intrinsic dust Small-Magellanic- 
Cloud-like (SMC-lik e) reddening in QSO spectra is discussed 
in some studies (e.g. lBaskin et aLll2013l ). For this single BAL 
QSO, the larger BAL-trough EW with flatter spectrum im¬ 
plies an important role for the ionizing continuum in the 
trough variation. 

The origin of BALs is commonly thought to be the re¬ 
sult of subtending a disk wind in the line of sight. The disk 
wind may be subject to different physical effects in different 
zones of its outflow velocity structure. It was found that the 
high-velocity components of the C iv BAL-trough change 
much more than the low-velocity part with the change of 
He II 1640 A EW, which is a measurement of the extreme 
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Lcont (xIO^ ergs s ' A ’) 

Figure 8. The EW versus the continuum luminosity for the C 
IV BAL (top), and for the Si iv BAL (bottom). The symbols are 
the same as in Fig. 7. 


UV continuum (Baskin et al. 2013). Here we divide the total 
trough EW into two parts, for high-velocity and low-velocity. 
For the C iv BAL-trough, the wavelength range is between 
1450 A and 1500 A (~ 20000 km s"^ - 10000 km s“^ ) for 
the high-velocity part, and between 1500 A and 1549 A (~ 
10000 km s“^ - 0 km s“^ ) for the low-velocity part. For the 
Si IV BAL-trough, the wavelength range is between 1340 A 
and 1355 A (~ 13000 km s“^ - 9700 km s“^ ) for the high- 
velocity part, and between 1355 A and 1400 A (~ 9700 km 
s“^ - 0 km s~^ ) for the low-velocity part. 

For the C iv BAL, the Spearman coefficients for the 
correlation between the divided BAL-trough EWs and the 
spectral index a are 0.64 (0.0042) for the high velocity part 
and 0.87 (0.000003) for the low-velocity part. The Spearman 
correlation coefficients are 0.75 (0.0003) for the high-velocity 
part of the Si iv BAL and 0.82 (0.00003) for the low-velocity 
part, where the values in brackets are the probabilities of the 
null hypothesis. It seems that the correlation between the 
spectral index and EW for the low-velocity part is slightly 
stronger than that for the high -velocity part. Tha t result 
is consistent with the study of (IZhang et al.l 120141 . where 
they found that the UV slope affects the low-velocity part 
of outflows much more. 


2 - r = -0.64 Pnull = 0.0038 
r = -0.47 Pnull = 0.05 


_ low velocity for C IV 

_ high velocity for C IV 











Lcont at 1300A 
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Lcont (xIO*^ ergs s~' A”') 
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Figure 9. C iv BAL EW vs the continuum, where the contin¬ 
uum luminosity is measured at 1300 A (top), at 1700 A (middle) 
and 2000 A (bottom), respectively. The black star is for the low- 
velocity part of the BAL-trough and the pink circle is for the 
high-velocity part of the BAL-trough. 


4.5 EW versus the continuum 

The dependence of the EW variability on QSO properties, 
such as the continuum luminosity, is still a question for de¬ 
bate (e.g., Gibson 2008; Filiz Ak et al. 2013; Baksin et al. 
2013). Fig. 8 displays the correlation between BAL-trough 
EW and the continuum luminosity for C iv (top panel) and 
Si IV (bottom panel). As in Fig. 7, blue squares denote the 
EWi data and red circles denote the EW 2 data. The correla¬ 
tion between the EW and the continuum is weak. The Spear¬ 
man coefficient is -0.19 (0.443) for C iv EW\, -0.21 (0.399) 
for C IV EW 2 , -0.24 (0.33) for Si iv EWi, and -0.22 (0.37) 
for EW 2 , where the values in brackets are the probabilities 
of the null hypothesis. The correlation between EW and 
the continuum is weak with a large probability of the null 
hypothesis. There is no significant correlation found in this 
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QSO. It is consistent with the study of iFiliz Ak et alJ (l2013l l 
who found no significant evidence for EW variability of C iv 
BAL-trough driven by the QSOs bolometric luminosity. We 
also calculated the Spearman coefficients between EW and 
the continuum for the high-velocity and low-velocity parts 
separately. The correlation coefficient between EW and the 
continuum is about -0.22 (0.39) (high-velocity) and -0.28 
(0.26) (low-velocity) for the C iv BAL-trough, and is -0.35 
(0.15) (high-velocity) and -0.22 (0.39) (low-velocity) for the 
Si IV BAL-trough, where the values in brackets are the prob¬ 
abilities of the null hypothesis. There is still no significant 
correlation between EW and the continuum for either the 
high-velocity or low-velocity parts. 

We also measured the continuum luminosity at 1300 A, 
1700 A, 2000 A, respectively (Fig. 9). The Spearman correla¬ 
tion coefficient between EW and the continuum luminosity 
(1300 A) for the C iv BAL is about -0.64 (0.0038) (low- 
velocity) and -0.47 (0.05) (high-velocity). For the continuum 
at 1700 A, the correlation coefficient is -0.366 (0.13) (low- 
velocity) and -0.3 (0.23) (high-velocity). For the continuum 
at 2000 A, the correlation coefficient is -0.059 (0.816) (low 
-velocity) and -0.063 (0.804) (high-velocity). The values in 
brackets are the probabilities of the null hypothesis. We find 
that the correlation becomes stronger when the continuum 
luminosity is measured in the bluer bands. The measured 
continuum luminosity in the bluer bands is more affected 
by dust extinction. This result is consistent with the the 
correlation between the BAL-trough EW and the spectral 
index. 


4.6 The maximum velocity of outflow 

The C IV maximum velocity of outflow Vmax for the 18 ob¬ 
servations changed by about 10% (Table 2). In Fig. 10, we 
show the relation between Vmax and the spectral index a. 
There is a moderately strong correlation between them with 
Spearman coefficient of 0.63 {Pnuii = 0.0005). This moder¬ 
ately strong correlation suggests larger extinction for larger 
maximum velocity of outflow. As we find a strong correla¬ 
tion between the BAL-trough EW and the spectral index, 
there should be a correlation between Vmax and the C iv 
BAL-trough EW, and indeed, the Spearman coefficient is 
0.78 {Pnuii = 0.00013). For flatter UV slope and the cor¬ 
respondingly larger BAL-tough, Vmax becomes larger. For 
the correlation between Vmax and the continuum(measured 
from ISOOA to 2400A), the Spearman coefficient is -0.33 
{PnuU ~ 0.174). A weak negative correlation seems to ex¬ 
ist between Vmax and the continuum with a probability of 
82.6%. The weak negative correlation between Vmax and the 
luminosity for this one QSO is not co nsistent with the pos i- 
tive correlation found by others (e.g. iLaor fc Brandd[2002lf . 
However, we should note that our correlation between Vmax 
and the luminosity is weak. More data are needed in the 
future for this study. 


5 SUMMARY 

The variability of the broad absorption lines is in¬ 
vestigated for a broad absorption line QSO, SDSS 
J022844.09+000217.0 (z = 2.719), with 18 SDSS/BOSS 
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Figure 10. The maximum velocity of outflow versus spectral in- 
dex a and the continuum (measured from ISOOA to 2400A). The 
Spearman correlation coefficient for the maximum velocity and 
spectral index a is 0.63 with PnuH = 0.0005. The Spearman cor¬ 
relation coefficient for the maximum velocity and the continuum 
is -0.33 with Pnull = 0.174. 


spectra covering 4128 days in the observational frame. The 
main conclusions can be summarized as follows: 

(1) The light curves for the BAL C iv EW, the contin¬ 
uum, and the spectral index are given. The EW for the C 
IV BAL-trough changes by about 25% (from 18.1 A to 22.9 
A). The a changes by about 260% (from -0.32 to -1.16). The 
continuum luminosity changes by about 38% (from 2.67 to 
3.69, in units of 10'*^erg s“^A“^). 

(2) Through the rms-to-mean ratio spectrum, the rela¬ 
tive flux change increases as the frequency increases, from 
about 10% on the red side, to about 15% on the blue side. 
The fractional variability is larger for the absorption trough 
than for the continuum and the emission lines. 

(3) There is a strong positive correlation between the 
C IV BAL-trough EW and UV spectral index a, as well as 
for the Si iv BAL-trough. The larger BAL-trough EW with 
flatter spectrum implies that the outflow has the effect of 
reddening the observed spectrum. The Spearman coefficient 
between the spectral index and BAL-trough EW for the low- 
velocity part is slightly larger than that for the high-velocity 
part. The strong correlation between the BAL-trough EW 
and the spectral index for this one QSO suggests that dust 
is intrinsic to outflows. 
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(4) There is no significant correlation between the C iv 
BAL-trough EW and the continuum in this QSO, nor for 
the Si IV BAL-trough. We do not hnd signihcant evidence 
for EW variability of the C iv BAL-trough driven by the 
QSO’s bolometric luminosity. We hnd that the correlation 
becomes stronger when the continuum luminosity is mea¬ 
sured in the bluer bands. The weak correlation between the 
BAL variability and the continuum luminosity for this one 
QSO implies that the BAL-trough variation is not domi¬ 
nated by photoionization. 

(5) There is a moderately strong correlation between the 
maximum velocity of the C iv BAL-trough and the spectral 
index a for this QSO. However, there is no signihcant cor¬ 
relation between the maximum velocity of outhow and the 
continuum. More observations are needed in the future for 
this study. 
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